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ABSTRACT 

Many Active Galactic Nuclei (AGN) are surrounded by gas which absorbs the radiation pro- 
duced by accretion onto the central black hole and obscures the nucleus from direct view. 
The dust component of the gas greatly enhances the effect of radiation pressure above that 
for Thomson scattering so that an AGN which is sub-Eddington for ionized gas in the usual 
sense can appear super-Eddington for cold dusty gas. The radiation-pressure enhancement 
factor depends on the AGN spectrum but ranges between unity and about 500, depending 
on the column density. It means that an AGN for which the absorption is long-lived should 
have a column density A^h > 5 x 10^'^ A cm~^ , where A is its Eddington fraction Lhoi/ L^dd, 
provided that iVn > 5 x 10^^ cm"'^ . We have compared the distribution of several samples 
of AGN - local, CDFS and Lockman Hole - with this expectation and find good agreement. 
We show that the limiting enhancement factor can explain the black hole mass - bulge mass 
relation and note that the effect of radiation pressure on dusty gas may be a key component in 
the feedback of momentum and energy from a central black hole to a galaxy. 
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1 INTRODUCTION 

Active Galactic Nuclei are powered by accretion onto a central 
massive black hole. The gas surrounding the nucleus, some of 
which provides the fuel for the accretion process, often obscures 
it from direct view. Indeed, the hard shape of the spectrum of the 
X-ray Back ground argues that most a ccretion onto galactic nuclei 
is obscured l lFabian & Iwasawalll999h . This is confirmed by deep 
Chandra and XMM-Newton imaging of the Sky with many AGN 
found to lie behind a column density of 10^^ — 1 0^"^ cm^^ or more 
jGiacconi et al ]|2002l : lBrandt & Hasingeill2005h . 

Much of the obscuring material must lie within the in- 
ner 100 pc, o r its total mass would be prohibitive (see 
iMaiolino & Risalitil2007l for a review). It is therefore part of the in- 
ner bulge of the host galaxy. Stars can also form from this gas, giv- 
ing a nuclear starburst. The gas is subject to the radiation pressure 
of the AGN, and can be ejected from the bulge if the nucleus be- 
comes too bright. Such AGN feedback may thereby remove the gas 
which fuels the nucleus and from which new stars form, so termi- 
nating the growth of both the central black hole and i ts host bulge. 
Simple calculations of when this occurs If abian 199?| jFabian et al.l 
l2002':'King'200l lMurrav et alj2005l : lFabian et al. 200ff) lead to the 
following relation between the mass of the black hole JV/bh and the 
velocity dispersion of the bulge a: 



good agreement with observations i Kormendv & Richston 1 19951: 



iMagorrian et al.lll998l : lGebhardt et al.ll2000l :l iFerrarese et al.ll200lh . 

The limit when radiation pressure ejects mass is an effective 
Eddington limit relying on absorption of radiation by dust, not elec- 
tron scattering. The radiation pressure is amplified or boosted by a 
factor A which is the ratio of the effective, frequency weighted, ab- 
sorption cross section for dusty gas, to that for electrons alone: 

0"T. 



A : 



(2) 



Mbh 
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Assuming a gas fraction / ~ 0.1, this gives an Mbh — o" relation in 



The dust absorption is greatest in the UV and the value of A de- 
pends on the spectrum of the AGN. The X-ray emission from the 
nucleus ensures that the gas and dust remain weakly ionized, and 
are effectively coupled by Coulomb forces so that pressure on the 
grains is shared with the surrounding gas . 

Boost factors computed for a standar d AGN spec t rum, using 
the radiation code CLOUDY are shown in iFabian et al. I ( 12006) and 
Fig. 1. They range from several hundred for low column densities 
and drop as the column density of gas increases until they approach 
unity when the gas becomes Compton thick, i.e. TVh ^ I/ctt ~ 
1.5 X 10^* cm~^ . Under the assumptions used, the main reduction 
of A with column density A'^h is merely due to the increased col- 
umn acting as a dead weight since the UV part of the spectrum is 
used up after traversing a column of a few 10^^ cm^^ . The outer 
'dead weight' is then pushed by the inner gas which experiences 
the full force of the radiation. 

In this paper we explore further the implications of the boost 
factor on the limit to the column density of dusty gas that can re- 
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log(N„/cm-') 

Figure 1. The radiation pressure boost factor A for dusty gas shown as a 
function of column density N^. The continuum sp ectrum is assumed to 
be tha t of high and low Eddington ratio objects (from lVasudevan & FabianI 
for the upper and lower lines, respectively. 
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Figure 2. Fig. 1 replotted in terms of column density versus 1/A. Long- 
lived absoiption from clouds near the centre of the galaxy should occur 
in the shaded region. Absorption from clouds and dust lanes can occur at 
higher values of A progressively further out. High column densities there 
would require prohibitive gas masses so we restrict such a region to A^h < 
5 X 10^^ cm~^ , marked by a horizontal Une. 



main close to an AGN without being ejected by radiation pressure. 
The key feature is to cast the inverse of the boost factor as the (clas- 
sical) Eddington ratio A: 

A-'^^ = A. (3) 

-t^Edd 

A cloud is long lived if A\ < 1. If the boost factor for a particular 
cloud is 100 then it will be blown away when the Eddington ratio 
exceeds 1/100. This determines a relation between the column den- 
sity of long-lived absorbing gas close to a nucleus and its Eddington 
ratiO; 

IPier & KrolilJ ( Il992l) have previously noted the role of the 
Eddington ratio when explaining the geometrical thickness of 
the torus around AGN by the action of radiation pressure. 
iHonig & BeckertI ( l2007b have estimated the effective Eddington 
limit for a torus around an AGN considering both smooth and 
clumpy dust distributions. 

Previous work on the luminosity dependence of absorp- 
tion has fo und that the inci d ence of absorp t ion d rops with 
luminosity (Uedaetal. 2003; La Fran ca et ay l2005l but see 
^welly & Page 2006, Treister & Urrv 200^. The 'receding torus' 
model jLawrencdl 199lUSimpsonll 1998t) is often invoked to explain 
why quasars may show less absorption. In that model the reduc- 
tion in absorbed objects is attributed to (heating) destruction of dust 
within the inner regions of a torus of absoring gas surrounding the 
nucleus. 

Here we consider the relation between luminosity and absorp- 
tion from the point of view of the effective Eddington limit on dusty 
gas. 



2 THE RELATION BETWEEN BOOST FACTOR AND 
EDDINGTON RATIO 

We determine the boost factor A as detailed in iFabian et alj ( |2006|) 
using CLOU DY and AGN spectral energy distributions obtained by 
I Vasudevari & Fabian (2007). A is obtained from absorption only 
(the input spectrum minus the transmitted one) and assumes that 
the gas is optically thin to the infrared radiation thereby produced. 
Trapping of radiation is assumed to be negligible and the ionization 
parameter ^ — Ljnr^ is arranged to be about 10 (thereby fixing the 
gas density for a given incident AGN flux). .Vasudevan & Fabiaiil 
(2007) find that the UV - X-ray spectral energy distributions 
(SEDs) of AGN depend on Eddington ratio, with much more ion- 
izing radiation - and therefore higher boost factors - occurring at 
higher A. We adopt mean SEDs for high (> 0.1) and low (< 0.1) 
A when computing A as a function of absorption column density. 
A^'h is plotted against A — A~^ in Fig. 2. The drift velocity of the 
grains relative to the gas is computed by CLOUDY and found to be 
low (about 0.1 km ^ where most of the radiative force is applied) 
justifying our assumption that the dust and gas are coupled. 

Long-lived absorbing clouds can survive against radiation 
pressure in the shaded region of the Figure. Clouds to the right 
of the dividing line, in the unshaded part, see the nucleus as above 
the effective Eddington luminosity and are thus ejected. Objects 
found in this region should be experiencing outflows and absorp- 
tion may be transient or variable. The Compton-thick objects never 
see the source as exceeding the Eddington luminosity and so can 
be long-lived at all Eddington ratios less than unity. Gas clouds in 
this regime are more likely to change because of star formation. 

Absorbed objects in the unshaded part at higher Eddington 
ratio should exhibit outflows, or the absorbing gas be far away from 
the nucleus where the retaining gravitational mass is much larger. 
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Th ey could for instance be associated with dust lanes, as envisaged 
Such absorption cannot be too large or the gas 
mass required would be prohibitive. We show, for example, a limit 
at iVn = 5 X 10^^ cm~^ in Fig. 2. 

Our prediction is therefore that absorbed objects should lie 
above the approximate dividing line given by A'h > 5 x 
lO^^Acm^^ , provided that TVh > 5 x 10^^ cm"^ . 



3 COMPARISON WITH DATA 

The simple model described above predicts that most highly- 
obscured AGN will be observed to have an intrinsic column density 
placing them in the shaded region of Fig. 2. 

In order to test our model, we have examined the absorption 
and black hole mass of several sam ples. The first is a com posite 
low redshift sample consisting of the lMarkwardt et ^ ( l2005h sam- 
ple of AGN detected by the Burst Alert Telescope (BAT) on the 
SWIFT satellite, and the Dadina (2007) sample of Seyfert nuclei 
observed by the BeppoSAX satellite. The former is an all-sky, hard- 
X-ray-flux limited sample in the energy range 14-195 keV where 
detection should be independent of column density, provided the 
source is not too Compton thick. The latter study presents an at- 
las of X-ray spectra for 39 Seyfert I and 42 Seyfert II nuclei over 
the energy range 2-100 keV. All objects in the combined sample 
have redshifts below 0.05. Black hole masses are estimated from 
veloci ty dispersions using the Af — u relation of iTremaine et al.l 
( l2002h . Velocity dispersions were obtained from the HyperLEDA 
online database. We obtained 2-lOkeV luminosities from the Tar- 
taruj3 database of ASCA observations, as it was not possible to 
extrapolate a 2-lOkeV lumin osity from the SWIFT luminosities 
provided in iMarkwardt et aL I {2005) without a value for the pho- 
ton index F. Finally, we convert to a bolometric luminosity us- 
ing the Eddington-rat i o-dep endent bolometric correction scheme of 
IVasudevan & FabianI ^20071) , using the black hole mass to estimate 
the Eddington ratio A, and assume that the emission is isotropic. 
We obtained velocity dispersions, and thus calculated masses, for 
23 objects in the composite sample. The results for this sample are 
shown in Fig. 3. 

In this Figure we see one object with a large column den- 
sity and high Eddington ratio in the unshaded region. This is NGC 
3783, and the abs orption is part of an outflowing warm absorber 
( lKaspietal.ll200ll) . The 3 objects with 22 > log iVn > 21 (just 
below the horizontal line) are, in order of increasing A, IC 4329A, 
which has an outflow and is s een almost edge on s o absorption 
may be from a distant dust lane jMarkowitz et al ."2006): N GC3516, 
which has variable absorption and an outflow ( Marko witz et al.l 
l2007h; and 3C120, w hich has a soft excess in XMM spectra 
( iBallantvne et al.l2004h so may have A'h overestimated in the value 
tabulated by Markwardt et al (2005) used in Fig. 3. 

We then used dee per samples. First we use the Chandra Deep 
Field South results of iTozzi et alj j2006) . These authors provide 
values for the column density A^h for each AGN together with 
intrinsic X-ray luminosities. We calculate black hole mass esti- 
mates using K-band mag nitudes from Szokoly et al. (2004) and the 
AIbh — Lk relation of Marconi & Hunll ( I2003i') . We impose a red- 
shift cut, requiring our objects to lie between redshifts 0.5 < z < 
1.0. We attempt to account for some evolution in the Mbh — Lk 
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Figure 3. Objects from the samples taken from the work of lMarkwardt et all 
(2005) and Dadina (2007) plotted on the A'h — A = iBol/^Edd plane. 
An estimate of the uncertainty in A = ieol/^Edd i** shown by the bar at 
top right. 



relation between that epoch and our own by incrementing the K- 
band magnitudes by unity before using the relation (i.e. account- 
ing for fainter bulge luminosities for the same central black hole 
mass in the past). This is broadly con sistent with the evolution ex- 
pected from Ivan der Wei et al J ( |2006|) . assuming negligible evolu- 
tion in the ratio Msteiiar/MBH. This produces estimates for A/bh 
in line with other measures. Using bolometric corrections as for the 
local sample, we estimate the Eddington ratios as before. This sam- 
ple yields 77 objects for which mass estimates could be obtained. 
The results are shown in Fig. 4 (top); objects with negligible col- 
umn density are marked nominally as having an upper limiting A'h 
of 10^^ cm^"^ (downward arrows). The range of Eddington ratios 
that we find fo r the CDFS sources is comparable to that found in- 
dependently bv lBabic et aLl ( l2007l) . 

Finally we used the Lockma n Hole results of iMainieri et al.l 

( I2OO2I) and iMateos et al] ( |2005|) . Mass estimates were calcu- 
lated using K-band magnitudes from Mai nieri et al. I ( l2002h , 
whereas X-ray lumi nosities and column densities were taken from 
iMateos et al ] l2005l) and bolometric corrections as before. The red- 
shift range used was the same as for the CDFS sample. There were 
13 objects satisfying these criteria, and the results are shown in 
Fig. 4 (bottom). 

The results from the deeper samples are qualitatively similar 
to the low redshift sample. For all three samples we find, as pre- 
dicted, that most objects lie within the shaded region. Note that 
since Chandra probes deeper than XMM, most of the CDFS ob- 
jects are at a lower value of A than those in the Lockman Hole. 

We show relatively conservative estimates for the uncertain- 
tiess involved to provide some indication of the validity of our re- 
sults. In the local Seyfert sample (Fig. 3), we use published errors 
on velocity dispersions and assume errors of 10 per cent on the Tar- 
tarus X-ray fluxes and also on the values of Nh provided. For the 
distant samples, we assume errors of 20 per cent on both K-band 
and 2-10 keV luminosities. In the case of the K-band, we note that 
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4 DISCUSSION 
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Figure 4. The A^h — ^ = ^Bol/^Edd p lane with objects from th e 
CDFS (top) using data from lXozzi etall i200d) and lSzokolv et alj j2004t)) , 
and the Lockman Hole (bottom) using data from lMainierieral 
and iMateos et al J j2005l) . There is one point missing from each plot with 
A < 10-4. 



Marconi and Hunt (2003) identify errors of ±0.1 on K-band mag- 
nitudes, translating to ~ 10 per cent in luminosity; we suggest that 
the 20 per cent value used here is therefore a reasonably 'safe' esti- 
mate of the error. For the X-ray band, one may expect variability to 
introduce significant uncertainties, but we suggest that 20 per cent 
error is again relatively conservative for most objects. Such error 
estimates lead to uncertainties of around ~ 40 percent in the Ed- 
dington ratios; however, this is small enough to clearly differentiate 
those lying within the shaded regions on Figs. 3 and 4 from those 
which are not. The typical error in A under these assumptions is 
shown in Figs. 3 and 4 by the bar in the top right. 



We see from the plots that most absorbed AGN with column den- 
sities A^H > W^^ cm"^ have an Eddington ratio A < 0.1, as ex- 
pected if radiation pressure on dusty gas is important. The observed 
objects are sub-Eddington in terms of radiation pressure on dusty 
gas and so the absorbing gas may be long-lived. The line defining 
the effective Eddington limit in the plots is determined for the mass 
of the black hole only. It shifts to the right (higher values of A) if 
larger masses associated with the stellar bulge are also involved. 
Since the lines drawn seem to apply to the samples, most of the 
absorbing gas must be at small radii (much less than 100 pc). The 
agreement between the expected location of absorbed AGN and the 
observed locations provides support for the idea that radiation pres- 
sure acts to blow gas away from galaxy bulges. If the radiation is 
isotropic then it can stem the growth of both inner bulge and black 
hole. 

The inner gas directly exp osed to the radiation p r essure 
may be unstable to c lumping ( iBIumenthal & Mathews! 1 19791 ; 
iHonig & BeckertI |2007|) , the details of which will not be pur- 
sued here. Most of the long-lived gas at column densities above 
10'^'^ cm"'^ is shielded by the inner gas and so need not be 
clumped. 

After blowing away the gas, AGN may decline in luminos- 
ity so creating unabsorbed sources at low Eddington ratios. There 
is also some uncertainty in the boundary between the various re- 
gions of the diagrams due to factors such as source variability and 
dumpiness in the absorption. The fact that most sources avoid the 
region of our diagrams above A^h ~ 5 x 10^^ cm"^ and to the 
right of the radiation pressure line demonstrates that variabilitity is 
not very important. 

Lower levels of absorption below 10^^ cm"^ can be long- 
lived at large radii in a galaxy, since the relevant gravitating mass 
there is due to the black hole and the bulge. The radiation limit in 
our Figures shifts to the right by a factor of Afbuigc/AfBH with 
all gas above the line bound to the bulge. The maximum boost 
A that can be obtained from radiation pressure on dusty gas is 
Cd/fT ~ 500 which means that the black hole is above the effec- 
tive Eddington limit for the whole bulge if A/Bn/A/buigc > 1/500. 

Consequently we envisage a scenario where a black hole 
smothered in gas could grow in a bulge in stages. It pushes the 
gas out to a distance in the bulge where the mass within that ra- 
dius is 500 times the black hole mass. (The boost factor increases 
as the column density decreases, so onc e gas starts to mov e out- 
ward it continues to do so, see Fig. 1 and lFabian et~ai]|2006l) After 
the accretion disk empties, the AGN switches off. If the bulge mass 
within the radius to which the gas was pushed exceeds 500 times 
the mass of the black hole, then the gas falls back in and the cycle 
repeats. Through accretion, the black hole mass and thus luminos- 
ity increases each cycle until is unable to retain the gas and it is 
pushed right out of the bulge. At this point 



MBH/Mbulgc - CTT/cTd ~ 1/500, 



(4) 



similar to the value found by Marconi & Hunt (2004) from corre- 
lating the observed properties of galaxies. 

Star formation from the gas in the galaxy during these cycles 
presumably leads to the bulge mass - velocity dispersion (Faber- 
Jackson) relation required such that equations (1), which acts lo- 
cally, and (2), which acts globally, agree. 

For much of the 'cycling scenario' envisaged above, the 
only acceptable range for bright unabsorbed objects would be at 
high Eddington ratios. From a large sample of AGN detected in 
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their A GN and Galaxy Evolution Survey (AGES). rKolImeier et alj 
( l2006h find A ^ 0.2. An important result from tlieir survey is tliat 
tliey should have been sensitive to unabsorbed AGN with lower 
Eddington ratios, but found none. This could in part be a selection 
effect due to absorption since the Chandra X-ray observations used 
are short, about 5 ks, which means that they are most sensitive to 
bright unabsorbed objects. The discussion of emission line strength 
for hi gh and low Eddington ratio AGN from lVasudevan & FabianI 
( 12007 ) could also be of particular relevance here, again implying 
that higher Eddington ratio AGN would be systematically favoured. 
The detected objects are in the unshaded part of our diagram, so 
have high A. 



5 SUMMARY 

Absorbed AGN are most commonly found at low Eddington ratios 
such that they are sub-Eddington for dusty gas. This agrees with 
the hypothesis that radiation pressure acting on dust is important in 
removing gas from galaxy bulges. In turn this leads to A/bh — cr 
and A/bh — Aifbuigc relations similar to those observed. 

Studies seeking to examine the evolution of absorbed AGN 
will need to include the dependence on Eddington ratio. 
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